Probing semi-macroscopic vacua by high fields of lasers 
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The invention of laser immediately enabled us to detect nonlinearities of photon interaction in 
matter, as manifested for example by Franken et al.'s detection of second harmonic generation and 
the excitation of the Brillouin forward scattering process. With the recent advancement in high 
power high energy laser and the examples of the nonhnearity study of laser-matter interaction by 
virtue of properly arranging laser and detectors, we envision the possibility of probing nonlinearities 
of photon interaction in vacuum over substantial spacetime scales compared with the microscopic 
scale provided by high energy accelerators. The hithertofore never detected Euler-Heisenberg non- 
linearities in quantum electrodynamics (QED) in vacuum should come within our reach of detection 
using intense laser fields. Also our method should put us in a position with a far greater sensitivity 
of probing possible light-mass fields that have been postulated. With the availability of a large 
number of coherent photons our suggested measurement methods include the phase sensitive (con- 
trast) imaging that avoids the pedestal noise and the scheme of second harmonic detection of photon 
nonlinearities in vacuum over a long co-propagating distance incurring resonance excitation. These 
methods carve out a substantial swath of new experimental parameter regimes of the exploration 
of photon nonlinearities in vacuum covering the force range from the electron mass scale to below 
neV. 

PACS numbers: 42.50.Pq, 42.65.Ky, 42.65.Es, 14.80.Va, 95.36.+X 



I. INTRODUCTION 



The invention of laser in 1960 constituted an introduc- 
tion of the possibility of coherent intense photon fields 
at the optical wavelengths. In fact immediately follow- 
ing it Franken et al.[l| observed the nonhnearity of the 
quartz crystal generated second harmonics of the laser 
at the 'high field' of 10^ V/cm. Many other nonlinear- 
ities of the media due to the high intensity of the laser 
have been discovered by subsequent years, as compiled 
in the late 1960's 01 ■ Above the field of manifestation 
of the nonhnearity of the second harmonic generation 
that Franken et al. encountered, neutral atoms may be 
directly ionized at the field called the Keldysh field [1], 
which is on the order of 10^ V/cm (i.e. an atomic scale 
of an eV potential over an A length), depending on the 
material. The material's nonlinearities arise from its po- 
larization under a large enough field of laser beyond the 
linear atomic field strength. See Figllja). Typically the 
restoring force of the electron to the rest of the atom 
saturates and can no longer match so much stronger ap- 
plied field. Relativistically strong lasers {i.e. the laser 
field is so strong to bring electrons to relativistic energies 
in an optical cycle) can induce the relativistic nonhnear- 
ity in electron dynamics of plasma, which is instrumen- 
tal in producing intense wakefields This is shown in 
Figlljb). Spurred by the promise and requirement of in- 
tense fields that can be sustained in a plasma wakefield 
excitation, since the invention of the CPA (chirped pulse 
amplification) [B*] the achievable laser intensity has been 
exponentially multiplying [6] . Along with this advance of 
most intense laser development comes a technique to fur- 
ther the available laser field such as in f?\, which in prin- 



ciple sees the possibility of reaching even the Schwinger 
field characterized by a quantum electrodynamic (QED) 
scale of an MeV (an electron-positron pair creation en- 
ergy) over the Compton length, i.e. 10^^ V/cm. Just as 
the material's ionization happens at the Keldysh field, 
the vacuum breakdown happens at the Schwinger field 
Q. This is depicted in Figlljc). Even before reaching 
this Schwinger field a host of nonlinear behaviors under 
intense fields are expected @. 

The study of matter and vacuum behaviors under in- 
tense fields these days is called high field science and 
reviews may be found in lol-fTH . In the following we 
suggest a class of possible investigations of vacuum in 
high field science, where we take advantage of the large 
amplitude of the (laser) electromagnetic fields and the 
relatively macroscopic scale of the field spatial scale (typ- 
ically micron) compared with the microscopic collider ex- 
periments (typically fm or less) . The most successful tool 
for exploring the fundamental nature in smallest struc- 
tures of matter has been that of collider, in which high 
energy (and thus high momentum) charged particles are 
produced and collide each other to probe the smallest 
scale with highest energies {Sx ~ h/p). The present ap- 
proach we suggest here may be contrasted to this high 
momentum approach, in that it explores semi-macro spa- 
tial scales and relatively lower energy processes with high 
amplitude (or fields) by 'exciting the constituent matter 
(or vacuum) and its structure\ 

We may wish to once again borrow the parallelism 
of the two alternative investigation methods of atomic 
physics by particle collisions and by laser nonlinear op- 
tics. The former explores the atom via Rutherford's 
method of the high energy beam streaming through and 
thereby scattered by the constituent matter deep in its 
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FIG. 1: Photon nonlinearities in media, (a) The nonlinearity 
arising from atomic anharmonic fields in high intensities; (b) 
The relativistic nonhnearity of plasma electrons in relativis- 
tically intense laser makes the harmonic oscillations of elec- 
trons (left) turn into anharmonic motions (right); (c) Even in 
vacuum photon nonlinearities occur when the field intensity 
approaches the Schwinger field, at which the intense field can 
polarize the vacuum anharmonically and eventually help vir- 
tual electrons and positrons turn into real particles, a clear 
manifestation of the breakdown of harmonic photon fields in 
the low intensity ordinary conditions. 



core. By this Rutherford discovered that the ordinary 
matter is composed of electrons that pervade the miit 
of the constituent matter, atom, while there exists a 
tiny core, the nucleus, that sharply scatters the incom- 
ing beam. The analysis of the scattering of the injected 
beam reveals the inner core property. Ever since this 
epoch making experiment, this approach paved the most 
successful way of the modern particle physics experimen- 
tal method, including the collider approach, that can ex- 
plore ever smaller spatial (and higher energy) structure 
of the constituent matter. On the other hand, the lat- 
ter optical approach available only after the invention of 
laser in 1960 was to excite the atomic electronic structure 
and carries out spectroscopy of the dynamics. The spec- 
troscopy reveals structure and its dynamical characters 
by palpating the atom by laser without pinpointed pene- 
tration into the core by the beam. This also includes the 
selective excitation of certain characteristic eigenmodes 
of the constituent matter structure. A good example may 
be the method called Laser Induced Breakdown Spec- 
troscopy (LIBS) When we study the vacuum itself, 
in stead of neutral atoms, this past experience guides us 
the distinction of the two approaches the high momen- 
tum vs. high field. As the main target of this research, 
in contrast to the collider's Rutherford approach with 
high momentum beam scattering, we are now introduc- 
ing an alternative approach of the method of 'exciting 
and probing the texture of vacuum' with high amplitude 
of photons. 



In order to probe such texture or polarization struc- 
tures of the excited vacuum, we are already treating vac- 
uum as a medium to be studied, rather than a given 
nothingness. For example, we posit that the phase veloc- 
ity of light is an experimentally meaningful measurable 
quantity. Intense photons fields (call it the pump laser 
or target laser) may provide nonlinearities even in vac- 
uum under a sufficient intensity (whose effects may be 
probed by another probe laser). Such was predicted by 
Schwinger [8] in the form of the generation of an electron- 
positron pair. However, manifestations of nonlinearities 
and thus the possibility of nonlinear spectroscopy of vac- 
uum should emerge even sufficiently below the Schwinger 
field threshold. (This may also correspond again to a sim- 
ilar phenomenon of the atomic ionization by intense laser 
at the Keldysh field and the emergence of atomic nonlin- 
earities far below this Keldysh field.) By this measure- 
ment we can investigate the dispersive and birefringence 
characteristics of photons in the excited vacuum. 

Let us consider the constituent's point of view, since 
the phase velocity shift can also be understood as a result 
of photon-photon interactions [13h17I] , as we understood 
the refractive index in matter through the photon- atom 
interaction at the level with atoms as the constituent. 
Photon-photon interactions in vacuum are quantum pro- 
cesses depending on the relevant frequency or mass scales 
of exchanged fields. Thus, the dispersive nature must 
be discussed based on the relevant frequency scale we 
introduce in experiments. As is known from particle 
physics, in 100 GeV scale we expect that the photon- 
photon interaction is caused by the heavy neutral bo- 
son exchange via higher order fermion loops to couple to 
photons |22| . As the frequency is lowered, lighter quark 
and electron loops may cause the photon-photon inter- 
action in the relevant quantum chromodynamics (QCD) 
and QED mass scales of lOOMeV and MeV. Below these 
scales there is no known mass scale relevant for photon- 
photon interactions. As long as we use lasers with the 
eV energy scale with the intensity below the Schwinger 
limit, therefore, the direct production of real constituent 
particles is not possible. However, it is possible to dis- 
cuss virtual vacuum polarizations. This was first intro- 
duced by the model through the Euler-Heisenberg effec- 
tive action [31] on the photon-photon interaction some 70 
years ago. Thus far, the real-photon-real-photon interac- 
tion has never been experimentally observed. Therefore, 
when we embark on the study of probing vacuum by op- 
tical laser fields, it is one of the first tasks to observe the 
QED interaction via the measurement of the phase ve- 
locity shift under the intense laser field due to this effect. 

In addition to exploring the nonlinear QED the cur- 
rent method may introduce a window through which we 
explore lighter mass scales well below MeV that have 
been speculated or hinted from particle physics [H, 
and cosmology These are related yet unobserved 

fields of axions and dark energy, for example. We note 
that the coupling of those light fields to matter must 
be extremely small. Otherwise, higher energy experi- 
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ments should have discovered these fields already, since 
light particles can be copiously produced. Here we recog- 
nize an advantage by introducing large amplitude rather 
than high momentum to search for these light fields that 
weakly couple to matter. The relevant mass scales of the 
exchanged fields between photons correspond to the force 
ranges. Therefore, the study of photon-photon interac- 
tions in the low frequency scale below optical frequency 
may unveil undiscovered semi-macroscopic forces as the 
basic building block enmeshed in vacuum. 

The focus of this paper is to discuss the experimental 
realization of a new type for the semi-macroscopic scope 
of vacuum by studying the property of nonlinearities of 
photons in vacuum. We shall discuss two examples of 
these. First we suggest an experimental technique of the 
phase contrast Fourier imaging to measure the phase ve- 
locity of photon under strong laser fields in section |TT] 
and mil This aims at a first verification of the nonlin- 
ear QED effect as pronounced by Euler-Heisenberg and 
beyond. Second we explore a realization of experiments 
to search for extremely light fields or long ranged force 
via the resonance interaction, employing quasi-parallel 
strong laser-laser interaction in section IIVI fVl and I VII In 
the conclusion section IVIII we summarize our approach 
with high fields of lasers based on a wider perspective 
by comparing our methods with those in particle collider 
physics and cosmology. 

II. PROBING NONLINEAR QED AND 
BEYOND 

Maxwell's equations suggest linear superposition of 
photons. On the other hand, the quantum electrodynam- 
ics (QED) indicates a weak interaction of photon-photon 
through the Feynmann's diagram of the box type, where 
a virtual electron-positron pair loops and four real pho- 
tons couple to the loop at the four vertexes as external 
lines (see FiglTl^b)). The photon- photon scattering cross 
section based on the box diagram is calculated by [28", '2^. 
The total elastic cross section in the center of mass sys- 
tem with the photon energy of uj is expressed as 

2 

where a = |- is the fine structure constant, rrip is electron 
he ' ^ 

mass and r, = a-^ - 2.8 x lO^^^ cm is the classical 
electron radius. For photons of ^ 1 eV, the cross 
section is Vi^^^ b. This is extremely small. The small- 
ness of this cross section arises from the electron-positron 
mass scale of the four propagators in the loop of the box 
diagram or from the short distance nature due to the 
corresponding Compton wavelength of electron-positron 
pair. We note that because of this smallness we have 
little 'noise', providing a pristine experimental environ- 
ment to search for something beyond QED. In reverse if 
we do detect any signals in photon-photon interactions. 



we are assured of something significant. When one is 
interested in observing the events of QED real-photon- 
real-photon scattering itself, the introduction of higher 
frequency photons as a probe beam onto a high intense 
optical laser target would have a greater probability due 
to the {JiLSf behavior than the optical-optical photon in- 
teraction. In such an asymmetric colliding system, even 
elastic collisions in the corresponding center-of-mass sys- 
tem (CMS) can be detected as frequency shifted scat- 
tered photons in the laboratory frame, depending on the 
collision geometry with reasonably high statistics. Fur- 
thermore, in such a setup, the non-perturbative nature of 
the intense field is expected to be important. For exam- 
ple, the catalysis of electron-positron pair production [s^l 
due to a higher tunneling probability from the Dirac sea 
may be tested, even below the Schwinger field Q. 

In what follows we discuss on the measurement of the 
phase velocity shift, where we focus on the optical-optical 
beam interaction. The strong electromagnetic field may 
modify the dispersion relation for photons. This effect 
is first discussed by Toll [l^ . If we detect the predicted 
velocity shift explained below, it amounts to the verifi- 
cation of the nonlinear QED effects in the perturbative 
regime. In the low frequency collision it is sufficient to 
describe the photon-photon interaction by the effective 
one-loop Lagrangian |3l|, HI] 

1 o? 

Li~ioop = ^^J^iPf^'^Pn' + 7{F,,Fn% (2) 

where Ff^^, — dAf^/dx" — dA^y/dx^ is the antisymmet- 
ric field strength tensor and its dual tensor F'^'^ = 
l/2e'"'^''Fxa with Levi-Civita symbol e'"'^'^ . It is well 
known that the forward scattering amplitude / with the 
dimension of length is related to the refractive index n 
via the Lorentz relation (ssj 

n(^) = l + -^7V/(w), (3) 

where N is the number density of scattering centers and 
/(w) is the forward scattering amplitude for light of en- 
ergy huj. The total cross section a is related to the ab- 
sorptive part of the forward scattering amplitude via the 
Optical Theorem, as follows 

47r 

a = — Im/(w). (4) 

UJ/C 

If we use a strong coherent electromagnetic field with 
large N, the real part of the forward scattering ampli- 
tude is expected to become large. This is because the 
coherent addition of the scattering amplitudes over large 
N is expected. Therefore, we are interested in the mea- 
surement of the refractive index shift under the high elec- 
tromagnetic field. It is not as effective to focus laser to 
cause non-forward scattering processes. The refractive 
index corresponds to the inverse of phase velocity. The 
derivation of these quantities in the linearly polarized 
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electromagnetic field target (so called crossed-field con- 
figuration where electric field E and magnetic field B are 
perpendicular with the same strength) is originally stud- 
ied in (T^ . [T5I and further derived from the generalized 
prescription based on the polarization tensor applicable 
to arbitrary external fields in [l^. This shows us 



v\\/c — 1 



v±/c = 1 — 



8 Zk 

14 Qf^ ^3 

45^ c5"fc2' 



(5) 



where v\\/c and v^/c are the phase velocities when the 
combination of linear polarizations of the probe and tar- 
get lasers is either parallel and perpendicular, respec- 
tively. The quantity mic^/H^ - 1.42 x lO^^J/m^ is the 
Compton energy density of an electron, k is the wave 
number of the probe electromagnetic field with the unit 
vector of k. The Lorentz invariant quantity Zk is defined 
as 



Zk = {k^F''-){kpF^), 



(6) 



and the relation to the energy density in the crossed 
field condition is 



Zk_ 
fc2 



e2(l + (fc.n))2, 



(7) 



with e — E — cB and fi = B x E. Thus the second 
terms in ([5]) show the deviation of the phase velocities 
of light W|| and v± are proportional to the field energy 
density normalized to the Compton energy density of an 
electron. 

The shift of the refractive index from that of the nor- 
mal vacuum is on the order of 10~^^ for the energy den- 
sity of IJ/fim!^. The refractive medium has the po- 
larization dependence, i.e. the birefringence nature. The 
difference in U|| and v±_ in ([5]) results from the first and 
second terms in the bracket of the effective one loop La- 
grangian in The UV limit (a; ^> cxd) of the disper- 
sion and birefringence under a constant electromagnetic 
field may be evaluated via the Kramers-Kronig disper- 
sion relation, as discussed in 17[. The phase velocity 



in both UV and IR limits is expected to be subluminal 
{vphase < c) under the QED field [H, The UV 

limit of the phase velocity is supposed to govern causal- 
ity. Therefore, it can be a fundamental test of a vari- 
ety of effective field theories in the IR limit whether the 
phase velocity in the UV limit extrapolated from that 
of IR is superluminal {vphase{oo) > c) or not. Thus far 
the dispersion relation from IR to UV is theoretically 
known only in the QED field [l^. However, there has 
never been data even in IR frequencies to date. It is 
important, therefore, for experiments to quantitatively 
verify or disprove the QED prediction. We note that 
the measurement of the refractive index in higher fre- 
quency may be sensitive to the part of the anomalous 
dispersion as discussed in [l^ and also the measure- 
ment of the electron-positron pair creation [l"9l - [2ll Isoj in 



strong electromagnetic field may be directly sensitive to 
the absorption or imaginary part. The Kramers-Kronig 
relation connects between the real and imaginary parts 
of the forward scattering amplitude or the refractive in- 
dex. Therefore, the systematic measurements of real and 
imaginary parts over wide frequency range may provide 
a test ground of QED and the Kramers-Kronig relation 
itself when it is applied to vacuum. 

Suppose then the detected dispersion and birefringence 
quantitatively deviate from the expectation of QED. This 
should indicate that undiscovered fields may be mediat- 
ing photons beyond QED. Scalar and pseudoscalar types 
of fields may appear via the first and second products 
in brackets of ([2]), respectively as we discuss section ITVl 
in detail. In addition if masses of scalar and/or pseu- 
doscalar fields are light, this long distance nature may 
enhance the coherence nature and thereby the amplitude 
of the forward scattering compared to the lowest order 
QED diagram. Therefore, the measurement of the bire- 
fringence and the comparison to the expectation from the 
nonlinear QED effects may be a general probe to inves- 
tigate the fundamental nature of vacuum. 



III. PHASE CONTRAST FOURIER IMAGING 
ON THE FOCAL PLANE 

We now consider an experiment where we create a 
high intensity spot by focusing a laser in vacuum and 
we probe its refractive index shift by a second laser. We 
call the first laser as the target laser, while the second 
as the probe laser hereon. The key issue is to detect an 
extremely small refractive index change resulting from 
the photon-photon interaction between the target and 
probe lasers. The conventional way that was already per- 
formed [131 and proposed [H, is based on the mea- 
surement of the ellipsoid caused by birefringence and that 
of the rotation angle of a linearly polarized probe laser 
by making it propagate for a long distance under a weak 
magnetic field [34j or electromagnetic field (s^UHl. This 
method has an advantage to enhance the phase shift by a 
long optical path without introducing costly strong tar- 
get electromagnetic fields. In the case of the constant 
magnetic field on the order of IT, one encounters the 
limit of physical detectability sensitive to the QED non- 
linear effect. In the case of strong electromagnetic field, 
one runs into the damage limit to store the strong tar- 
get field within a cavity over a long time. On the other 
hand, if we can directly utilize the local nature of vacuum 
by tightly focusing an intense laser pulse and measuring 
the lensing effect of vacuum on the pulse-by-pulse basis, 
there is no physical limit in increasing the intensity of 
the laser pulse until vacuum itself breaks down. In order 
to increase the shift in the refractive index corresponding 
to the inverse of phase velocities in ([5]), «.e. the intensity 
of the target laser field as expected through ^ and ([7]) , 
it is necessary to use a focused laser pulse by confining 
a large laser energy into a small spacetime volume. This 
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FIG. 2: Diffraction images at the far distance from the thin 
wire target when the Gaussian laser is shot on the wire. The 
left figure is the case there is no wire. The right figure is the 
case when the thin wire target is horizontally arranged. 



causes a locally different refractive index along the tra- 
jectory of the target laser pulse in vacuum. A variation 
of the refractive index arises over the high intense part 
and the rest of vacuum. If the probe laser penetrates into 
both parts simultaneously, the corresponding phase con- 
trast should be embedded in the transverse profile of the 
same probe laser. Our proposal is to directly measure the 
phase contrast and to determine the absolute value of the 
refractive index change by controlling the combination of 
polarizations of the probe and target laser pulses. This 
should result in the birefringence as expected in ([5]). 

We need to detect the extremely small shift of phase 
velocity by the target-probe interaction. For this we also 
need an intense probe laser in order to enhance visibil- 
ity. However, if we utilize the conventional interferome- 
ter techniques providing a homogeneous phase contrast 
over the probe laser profile, such small refractive index 
changes are hard to detect. This is because the result- 
ing intensity modulation is always on top of the huge 
pedestal intensity and the contrast of the modulation to 
the pedestal is extremely small. Any photo device cannot 
be sensitive to a small number of photons spatially mod- 
ulating under the pedestal intensity beyond IJ (~ 10^^ 
visible photons) due to the limit of the dynamic range on 
the photon intensity measurable by a camera pixel with- 
out causing saturation on the intensity measurement. On 
the other hand, broadening the dynamic range by lower- 
ing the gain of the electric amplification of photo-electron 
degrades the sensitivity to the small number of the spa- 
tially modulating photons or the sensitivity to the small 
phase shift. Therefore, we need to invent a method that 
can spatially separate the weakly modulating character- 
istic intensity pattern from the strong pedestal. 

In order to overcome this difficulty, we suggest uti- 
lizing the inhomogeneous phase contrast Fourier imag- 



ing on the focal plane by focusing the probe laser via a 
conceptual lens. In this the physically embedded phase 
contrast on the transverse profile of the probe laser am- 
plitude is Fourier transformed on the focal plane. The 
intensity pattern on the focal plane has a preferable fea- 
ture that the characteristic phase boundary makes the 
intensity profile expand outer regions far from the fo- 
cal point, whereas a Gaussian laser with a homogeneous 
phase converges into a small focal point with the smallest 
beam waist. It is instructive to illustrate the character- 
istic nature of the diffraction pattern from a wire-like 
target shape as shown in Figl^l Here the intensity pat- 
tern at a far distance known as the Fraunhofer diffraction 
limit is shown in the case when a Gaussian laser beam 
is shot on the thin wire target. The intensity pattern at 
a far distance can be understood as the Fourier trans- 
form of the wire shape approximated as a rectangular of 
2pL X 2v. It is well known that a lens component by it- 
self also has the equivalent effect to produce a diffraction 
pattern corresponding to the exact Fourier transformed 
image of the original shape of the refractive medium on 
the focal plane at a finite distance (see [s^ for instance) . 
In order to understand the diffraction image, we may re- 
fer to Babinet's principle which states that the diffraction 
pattern from an opaque wire is identical to that from a 
slit of the same size and shape. The Fourier transform of 
such a rectangular slit is expressed as 

/ siii(^iw^\ ^ / sin(z^\ ^ 

V IJLLOx ) \ VUJy J 

where ujx — and ojy ~ are the spatial frequencies 
for the given position (x, y) on the focal plane with the fo- 
cal length / and the wavelength A. In the case of the slit 
with ^ ^ the rectangular profile on the focal plane be- 
comes orthogonally rotated thinner line shape with oscil- 
lations included in the line (see Figl2] (right)). This is be- 
cause the narrower the slit size is, the smaller the spatial 
frequency in that direction becomes. On the other hand, 
the shape of Gaussian distribution without the wire or 
slit is unchanged, because the Fourier transform of the 
Gaussian distribution corresponds to the Gaussian dis- 
tribution (see FigUl (left)). This is the key feature that 
drastically improves the detectability to the small phase 
shift by sampling only outer parts far from the Gaussian 
part. This may also be interpreted as the counter-concept 
to the conventional spatial filter, where outer parts are 
eliminated to maintain a smooth phase on the transverse 
profile of the Gaussian distribution. 

Given the intuitive picture above, a quantitative for- 
mulation of our proposed method is presented as follows. 
First, let us define the geometry of the laser intersection, 
as shown in FigI3] (left), where the tightly focused target 
pulse with time duration tj propagates along the Z-axis 
and the probe pulse with the larger profile and longer 
time duration Tp propagates along the z-axis tilted by 
6 from the Z-axis. In this figure the wavefronts of the 
probe pulse are drawn successively with time step tj un- 
der the condition Tp > tj. Since the constant phase shift 
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FIG. 3: Geometry of the embedded phase contrast in the 
probe pulse. 



is embedded only during Tt, the optical length with the 
constant refractive index shift along the z-axis perpen- 
dicular to the wavefront is limited to ctj. This is inde- 
pendent of 0. In this geometry after the penetration of 
the probe laser pulse, the profile of the probe laser on the 
X — y plane contains a trajectory with the constant phase 
shift 6 along the projection of the path of the target laser 
on the probe wavefront, as shown in Fig |3] (right). 

In order to discuss the amount of the phase shift, we 
need a concrete geometry at the diffraction limit of both 
target and probe lasers. Let us briefly review the laser 
profile at the diffraction limit. The Gaussian profile is 
a basic constraint in typical laser fields, where the aper- 
ture of a lasing material has a finite size in the transverse 
area. The solution of the electromagnetic field propaga- 
tion in vacuum with the Gaussian profile on the trans- 
verse plane with respect to the propagation direction z is 
well-known j40| . The electric field component in spatial 
coordinates (a;, y, z) is expressed as 



E{x,y,z) oc 

Wo 



w{z) 



exp ■ 



'i[kz — ri{z)] — 



ik 



wizY 2i?(z) 



,(9) 



where k — 2tt / \^ r — yx^ 



y^ , wa is the minimum waist 



which cannot be smaller than A due to the diffraction 
limit, and other definitions are as follows: 



w{z) 



Wo 



zr 



R = 



•qiz) = tan M 1 + — 

ZR 



(10) 



(11) 



(12) 



ZR 



A 



(13) 



Let us express the phase shift 5 in the limit z = cTf < 
zjj, where we assume an almost flat wavefront as indi- 
cated by pO)) and dTTI) in the vicinity of the diffraction 
limit 



5 = —SnqedCTtLpt{Xp,yp), 
Ap 



(14) 



where the subscripts p and t denote the probe and target 
quantities, respectively, Sn is the refractive index shift, 
CTt is the path length with effectively constant phase shift 
S over the crossing time Tt and ipt{xp,yp) is a weight 
function to reflect the path length difference depending 
on the incident position with respect to the target proflle 
expressed as a function of position {xp, yp) on the trans- 
verse plane of the probe laser. Based on ([5]), ((H]) and ([7]) 
we parametrize the refractive index shift as 



5nqed = C^o(l - cos^)- 



(15) 



where C is 4 and 7 for the polarization combinations || 
and _L, respectively, No is the coefficient to convert from 
energy density to refractive index shift defined as iVo = 
A^g. ^ 1.67 X 10^12 [^fm^/J], is the incident angle of 
the probe pulse with respect to the propagation direction 
of the target pulse as depicted in Figl^l Et is the energy 
of the target pulse in [J], and nwof-CTt is the volume in 
[/im'^] for the given target profile with the minimum waist 
from (|10p . By substituting (ITSI) into we obtain 
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the experimentally convenient expression for S 



An 



^)—^Vt{xp,yp). 



nw, 



(16) 
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We are interested in an application in the limit wot < 
wop for the width and en < wop for the depth of the 
embedded phase shape. We then take the approximation 
ipt{xp,yp) ~ 1 to simplify the following argument (if nec- 
essary, we may restore the target profile ft{xp, yp) based 
on the precise profile of the target laser refiecting actual 
experimental setups). In this limit we approximate that 
the target profile has a rectangular shape with the size of 
2fi X 21^ inside which a constant phase shift is assigned, 
where the effective slit sizes are defined by the transverse 
sizes of focused laser beams through the relation fi ^ wop 
and p wof We note that the optical length cTt with 
phase shift S of the probe laser is eventually canceled 
out based on This indicates that we have many 

choices on Tt as long as the conditions CTt < z^t ^^nd 
Tt < Tp ^ ziit/c are satisfied. 

We introduce the window functions rec and rec as fol- 
lows: 



rec{fi, v) : 
rec{^, v) 



1 for < ^ n \y\ < u 

for |x| > ^ U \y\ > v 

for \x\<pr\ \y\ < V 

1 for > /X U \y\ > V 



(17) 
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This window provides a unit region of a constant phase 
which is employed for arbitrary phase maps composed of 
a collection of the unit window cells. If we determine the 
phase on the pixel-by-pixel basis for a given camera de- 
vice, the rectangular shape as a minimum unit cell should 
be a natural choice. 

The intensity distribution at the focal plane is deter- 
mined as a function of the peak amplitude of the probe 
pulse Aop, the wavelength Xp and the focal length fp for 

a given Gaussian beam profile of Aope~'^^^«^'^°'^ as the 
probe laser. The linearly synthesized amplitude at z af- 
ter crossing with the target laser pulse is expressed as 

ipixo, yo) = arec{n, zy)e"°'^''o+^o) -|- l3rec{p, i/)e"°^^"+^o\ 

(18) 

where a and /3 are the plane waves at the point z after 
the penetration of the probe laser. The function a with 
the local phase shift 5 caused by the local refractive index 
shift and (3 without phase shift are defined as 

a = Aope''-'''+'\ 

/3 = Aope^'"'. (19) 

The Fourier transform F of the synthesized amplitude 4' 
on the focal plane (x, y) at z is expressed as 

F{i^{xo,yQ)} = 
aF{rec(/i, iy)e"°("^o+!/o)} + j3F{rec{fi, t/)e-°(^o+yo)} 

= (a~f3) r I dxodyoe~''^''°+y»'>e-''^'^'''="+'^yy°'> + 

/oo poo 
/ dxodyoe-''^'=»+y''>e''^'^y'=''+'^yy'>'>, (20) 
-ocJ —oo 

where we define (uJx,^v) = 2;, -^^y). We introduce 
coefficient Csig for the first term in (pO)) containing the 
information on how much the phase shift, namely signal, 
is localized resulted in the photon-photon interaction de- 
fined as 

Csig i^x, ^y) = 

dxQe''''-'^ cosily xXq) I dy^e^^y^ cos{ujyya). (21) 

- p, J —u 

We also define coefficient C^kg for the second term of (PO)) 
which corresponds to the background pedestal Gaussian 
part as 

TT (^x+'^y) 

Cbkg{^x,^y) = -e ■ (22) 

Therefore, the Fourier transform becomes 

F{V'(xo,2/o)} = (a - (i)Csig{uJx,i^y) 

+l3Cbkg{i^x,^y)- (23) 

By substituting (dH), (ED) and ^ into the inten- 
sity pattern at the focal plane is expressed as 

\tl;{ujx,ujy)\'^ = 
(^)2{2a.,(C..g - CbkgKl - cos 5) + Clg}. (24) 



Equation (l24l) indicates that this method eventually cor- 
responds to an interferometer via the cross term of 
2Csig{Csig — Cbfcg)(l — COS (5). This interferometer dif- 
fers from the conventional one in that the modulating 
part due to the phase shift d can be spatially separated 
from the confined strong Gaussian part C'^^.g due to the 
nature of Csig part. Therefore, in principle, our method 
of sampling only peripheral intensity modulations caused 
by non-zero 6, provides a high signal-to-pedestal ratio cir- 
cumventing the most intense focal spot. This is demon- 
strated in llllBI See Figl5]and|6l 



We note that this method is found to be similar to but 
distinct from the idea in [37[, where two intense target 
laser pulses are treated as a matterless double slit and 
the interference between spherical waves from these slits 
is discussed as a signature of the photon-photon interac- 
tion. In their proposal the physical diffraction is caused 
by laser-laser interaction itself. Our method rather lets 
the target laser cause the refractive phase shift in the 
probe laser as indicated in Fig|4l This phase shift is em- 
bedded in a refracted nearly plane wave in the forward 
direction of the probe laser as explicitly formulated in 
([TS]) and (ITOl) . We then set a lens to the right of inter- 
action between the target and probe lasers as shown in 
Fig HI The diffraction or Fourier transform in our method 
is incurred by the added phase of the lens component and 
the spherical wave propagation from the lens to the fo- 
cal plane. The advantage of our method is the enhanced 
detectability of small phase shift on the pulse-by-pulse 
basis, as is demonstrated in section ITlIBI due to more ef- 
ficient collection of photons by the lens effect with the 
simpler target geometry. On the other hand, a disad- 
vantage is the deviation from the ideal lens phase after 
embedding the physical phase on the probe laser. How 
to correct this kind of background phase fiuctuations is 
discussed in the following subsections. 



We also note that the Fourier image on the focal plane 
provides only the absolute value of phase shift 5. Because 
of 1 - cos (5 ^ with (5 < 1 in the proposed 

method has no sensitivity to the sign of phase shift in 
the case of almost homogeneous background phase. On 
the other hand, if a special offset phase tt/2 is locally 
added along the path of the target laser before the final 
focusing to have the Fourier image, such setup may re- 
cover the sensitivity to the sign of phase shift, because 
1 — cos(^ + tt/2) in becomes ^ 1 + S. From physical 
point of view, it is important to discuss whether phase 
shift is positive or negative, since it refiects the dynamics 
of local interaction. Furthermore, this has a definite ad- 
vantage to enhance the signal due to the proportionality 
to 6 compared to 5^ for the extremely small 5. There- 
fore, the implementation of such local offset phase on the 
probe laser in advance should be a part of experimental 
task. 
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A. Template analysis for local phase reconstruction 

In actual experimental situations it is unavoidable to 
contain local phase fluctuations inside the probe pulse on 
the pulse-by-pulse basis even in the absence of the phys- 
ical signal. Deviations from the ideal Gaussian mode [39] 
and local phase fluctuations due to optical elements in 
the path of the probe laser may be sources of local phase 
fluctuations. Compared to the physical phase shift 6 from 
the nonlinear QED effect, these fluctuations are expected 
to be much larger. However, if the local phase map 4>{X) 
on each probe pulse is a priori measured as a function of 
position X = (xq, yo) on the transverse plane of the probe 
pulse, we may be able to correct the effect from the back- 
ground fluctuations. In the next subsection IIII Bl we dis- 
cuss how to measure the discrete phase map (jfi = 4>{Xi) 
on the pulse-by-pulse basis in detail, where i denotes a 
discrete position on the transverse plane at which 6 is 
embedded. Here we focus on how to determine phase 
shift S from the combination of experimental data on the 
phase maps under such background fluctuations on the 
pulse-by-pulse basis. 

Let us generalize the discussion of ([20]) . The integral 
range in (IT71) included in the first of ^U\\ can take any 
shape and size in general. We replace the rectangular 
region rec with the region Ri = R{Xi) where a con- 
stant phase is mapped within the region for the given 
discrete position Xi. By denoting the spatial frequency 
as = {ujx,^y) = {2Trx /{fpXp),2ny/{fpXp)) for the po- 
sition on the focal plane (x, y) with the integral kernel 
f{W,X) = e-"(=^o+yo)e-^('^-^t)+"!'2'o), the Fourier trans- 
form including the local phase fluctuations 0i is expressed 
as 

i^iW;^) = F{ij{X;4>)} 

Nx „ 

= V{a(</),)-/3} / dXfiX,W) 

Nx 

= ^{a(</),)-/3}/.(M^)+/3/ooW, (25) 

i 

where Nx is the number of regions on the transverse 
plane at z, a{(j},) = Aoe-*('''^+'^'), /3 = ^oe"*'=^ h{W) = 
Jj^^dXf{X,W), and Io.{W) = JZ,dXf{X,W). We 
note that this expression corresponds to the regional cut 
and paste of a Gaussian laser amplitude; i.e., cutting a 
region with a phase determined from (3 at z and paste 
the same region by adding in a{(l>i). 

Given (/'(X^) on the pulse-by-pulse basis, we can 
numerically calculate the real and imaginary parts of 
^p{W;4>) as follows; 

Nx 

Reij{W; (j)) ^ ^{cos(fcz + </>,;) - cos{kz)}I^{W) 

i 

+ coa{kz)Iao{W), 



Nx 

lyji4j{W; (/)) = ^{sin(fcz + 0,) - sm{kz)}I,{W) 

i 

+ sm{kz)Ioo{W). (26) 

The estimated background intensity pattern Ibkg{4') on 
the focal plane with the phase fluctuations without 
physical phase is given by 

hkgiW; 0) = {Re^iW; (b)}^ + {Im^(VF; (27) 

We now embed a template of the local physical phase 
shift Si = 6{Xi) as well. The phase shift 5i can be eval- 
uated from geometry of the energy density profile, i.e. 
the refractive index distribution inside the intense tar- 
get laser field. This can be fixed from the experimental 
design on the focal spot a priori. We can monitor if 
the center of the spot is surely stable and further cor- 
rect its deviation from the fixed geometry of the target 
laser. Given Si, we have to only replace the phase by 
<j)i ^ (pi + uSi with a constant parameter k as follows 

Ibkg+sig{W;(j) + kS) = 

{ReV^W^; (j) + K6)y + {lin^(W; (j) + k5)Y, (28) 

where bkg -\- sig refers to the fact that the physical phase 
is embedded in the background phase. 

Given the measured intensity pattern Imeas on the fo- 
cal plane per probe pulse focusing, we define with (j28p 
as a function of k 

^■meas {Wj )+ hkg+sig {Wj ;(!) + kS) 

where Nyy is the number of sampling points on the fo- 
cal plane and j runs through all discrete positions on 
the focal plane except for the most intense region around 
the focal point. We then determine the parameter k by 
minimizing the on the pulse-by-pulse basis within the 
acceptable accuracy. The reason why we have to deter- 
mine K on the Fourier image on the focal plane is due 
to the experimental constraint that we cannot sample 
the most intense part at the focal point. Because of the 
information loss at the focal point corresponding to the 
lower spatial frequency part, we cannot reconstruct the 
original phases of the probe laser before focusing by the 
inverse Fourier transform for the lens effect. 

We note that determining k corresponds to the mea- 
surement of the absolute value of phase shift S. As we 
have discussed in the last paragraph of section IIIIl if the 
offset phase 7r/2 is locally added to the known physical 
template before the focus by the lens, it provides the sen- 
sitivity to the sign of phase shift. This also has the advan- 
tage that the intensity change of the physical signal on 
the characteristic Fourier image is drastically enhanced 
due to the proportionality to i5 compared to 6^. Such 
local phase may be implemented on the probe laser by 



9 




FIG. 4: Conceptual experimental setup for the suggested 
phase contrast Fourier imaging. At the crossing point between 
the probe and target lasers, the target laser causes the shift in 
the index of refraction, which amounts to the refractive phase 
shift S embedded in the probe laser as explained in Fig[3] 
After the expansion from the crossing point, this probe laser 
goes through the lens and the signal path on the right. The 
diffraction or Fourier transform is incurred by the added phase 
of the lens component and the spherical wave propagation 
from the lens to the focal plane as demonstrated in Fig[2]and 
later in Fig[5]and[6l 

utilizing the photorefractive crystal. The photorefractive 
effect [40| causes a static local refractive index change by 
supplying an external electric field onto the crystal. If we 
embed a rectangular shape with different refractive index 
on the surface of the photorefractive crystal by mimick- 
ing the shape of the target laser trajectory in vacuum, we 
may be able to implement the local phase offset 7r/2 onto 
the physical template in advance by making the probe 
laser propagate through the crystal before the final fo- 
cusing by the lens. If this is the case, we can use a much 
weaker probe laser. In following subsections, however, we 
consider the case where such special offset phase is not 
implemented; instead, a sufficiently intense probe laser 
as well as the target laser is available in an experiment. 



B. Suggested setup and the Fourier image 

Figure U illustrates the conceptual experimental setup 
for the phase contrast Fourier imaging on the focal plane. 
The setup consists of two parts; the signal path (SP) and 
the calibration path (CP). The part SP is to perform the 
phase contrast Fourier imaging to measure a physically 
embedded phase shift by the probe-target laser interac- 
tion. Combined with SP, CP is to provide phase maps 
of static optical components (j)opt{Xi) and the phase map 
(j)pis(Xi) of probe laser brought on by instabilities in the 
upstream laser systems on the pulse-by-pulse basis, as we 
have discussed in llll Al 



Within SP, both the target and probe laser beams are 
focused with different waist sizes Wto and Wp^^ with the 
i^-number. Ft — ft/d and Fp — fp/d for the common in- 
cident beam diameter d, respectively. The lens size must 
be larger than d. They are crossed each other at the mini- 
mum waist point where the wave fronts in both beams are 
close to flat in the case of the ideal Gaussian profile beam 
with i? = oo in pT|) . We assume that the target waist 
Wot is smaller than the probe waist wop which embeds 
a phase contrast within the amplitude on the transverse 
profile of the probe laser. 

The crossing angle 9 has no impact on the optical 
length with the shift of the refractive index. However, 
it is important to control the amount of the refractive in- 
dex shift, depending on the relevant dynamics for photon- 
photon interactions. In the case of QED the setting 6* ~ tt 
has the advantage to increase the photon-photon center 
of mass energy, as discussed in We design the tar- 

get geometry such that the Rayleigh length of the target 
pulse Zfit is approximated as fi in (I17[) and the target 
beam waist Wto is approximated as v. Accordingly the 
beam waist of the probe pulse Wp^ should coincide with 
Zflj within the probe duration time Tp = zj^t/c. 

After embedding the physical phase shift onto the 
probe pulse at the minimum waists of both beams, the 
probe laser is expanded by expansion factor Re to get 
the probe beam waist oi Wp^. The expanded probe laser 
is refocused with the focal length fp^. It is important to 
notice that the amount of the embedded phase contrast 
is independent of the rate of magnification of the trans- 
verse profile of the probe laser. This implies that we have 
degrees of freedom to design the final focus to produce a 
Fourier image with the sufficient number of photons per 
camera pixel that enables the pulse-by-pulse analysis, as 
we see below. 

The procedure to determine k by the physical inter- 
action is as follows. First, we determine constant phase 
biases (f>opt{Xi) caused by all optical elements included 
in both SP and CP by performing the phase contrast 
Fourier imaging at the ends of CP and SP. We may as- 
sume that (poptiXi) is stable over many pulse injections. 
The phase (j>opt{Xi) can be measured with a CW laser 
over a long period in between target-probe beam cross- 
ings during an experiment. Since a less intense single 
mode Gaussian CW laser can provide a homogeneous sta- 
ble phase as the average value, we can, in principle, de- 
termine 4>opt(,Xi) with high accuracy based on the phase 
contrast Fourier imaging integrated over the long time 
period at CP and SP by sharing the common CW laser. 
The template analysis discussed in IIII Al can also be ap- 
plied to this purpose. By assigning a square shape to 
the region Ri representing a pixel instead of the signal 
in (P^. we estimate k for each Xi. The number of pho- 
tons at a point Wi on the focal plane contains convoluted 
phase information on the amplitude from all points on 
the transverse plane of the probe Nx as seen from (|25l) . 
Therefore, as long as the number of sampling points on 
the focal plane Nw is larger than that on the transverse 
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probe profile Nx^ we can, in principle, determine a set 
of 4>opt{Xi) from by scanning k for each Xi over 
the expected range of the phase variation. The achiev- 
able resolution of the phase reconstruction depends on 
the scanning step on k, in the test. As we see later in 
TabUl FigEl and FigEl the phase contrast Fourier imag- 
ing at a focal plane achieves the sensitivity to the order 
of ~ 10~^ in the physical phase shift by sampling only 
outer regions on the focal plane with a several-photon- 
sensitive photo device. Therefore, we can introduce the 
same resolution step to determine n for 4>opt[Xi) as the 
analysis of the physical phase shift 5{Xi). We may mea- 
sure the initial coarse phase map <j)opt{Xi) a priori by 
a typical phase resolution ^ A/ 100 with a commercially 
available wave front sensor. In such a case we need to 
repeat the two dimensional Fourier transform about 10^ 
times for scanning k for each Xi to reach the same phase 
resolution as ^ 10^^ starting from the initial coarse res- 
olution. Accordingly, a proper computing power is nec- 
essary to solve the local phases over a set of (j)opt{Xi). 
Second, we can also apply the same procedure to deter- 
mine a set of 4)pis{Xi) on the pulse- by-pulse basis after 
mapping the static phase (popti^i) without the physical 
template 5{Xi) within CP. Finally, we create a phase map 
by adding (t)opt{Xi) and (j)pis{Xi) within SP. We then de- 
termine K caused by the physical template S(Xi) based 
on with the same phase resolution step when we 
determined (poptiXi) and (ppuiXi). 

As a demonstration, FiglSftop) illustrates the intensity 
distribution on the focal plane in SP with an arbitrary 
unit on the contour height in log-scale by sampling val- 
ues with 50/im steps in both x and y-axes. In this figure, 
only the physical template 5{Xi) is embedded. The pa- 
rameters used for Figl5]are summarized in TablH where 
the parameters for both target and probe lasers, the em- 
bedded physical phase shift due to the nonlinear QED 
effect, and focusing parameter after the beam expansion 
to obtain the Fourier transformed intensity distribution 
in SP are specified. Figure [5] (bottom) shows the same 
distribution without the physical phase shift in the 3- 
D plot, where the pedestal Gaussian distribution is well 
confined within an 20 x 2Q^Ta^ area. 

Figure [6] (top) and (bottom) show the integrated num- 
ber of photons in a single target-probe crossing per 50/im 
X 50/im pixel along the y and x-axes on the focal plane, 
respectively. These correspond to the photon yield per 
camera pixel in single pixel line along the y and x-axes, 
respectively in Fig IS] The plotted range is extended over 
±1 cm. In each of FigjBl the number of photons in 
the peak region confined in a pixel without phase shift 
is overlaid. This indicates that the Gaussian pedestal 
part is well confined and the modulation from the phys- 
ical phase shift is expanded to the outer directions on 
the focal plane. Therefore, by sampling only the pe- 
ripheral surrounding the focal point on the focal plane 
without touching the most intense pedestal, we should 
be able to count the sufficient number of photons per 
pixel on the pulse-by-pulse basis. We may introduce a 



With phase shift 




x[nm] 



Without phase shift 




-50 .50 



FIG. 5: Intensity profile in unit of the number of photons 
per single shot focusing on the focal plane by sampling the 
values with 50/im steps in both x and y axises. The relevant 
laser parameters are summarized in Tab|T] The top figure is 
the intensity value per point with physical phase shift based 
on (|24|> . The color map is normalized to 10'' by truncating 
the intensity at the focal point to zoom the signal part in the 
color map. The bottom plot shows the same quantity with- 
out physical phase shift in 3-D to show how much the ideal 
pedestal Gaussian intensity profile is localized quantitatively 
in the x — y plane without the truncation of the intensity at 
the focal point. 
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Target laser parameters 


Probe laser parameters 


Tt = 10 fs 


Tp = ZRt/c = 12 fs 


Et = lOkJ ( 4.03 X 10^^ photons ) Ep 


= lOkJ ( 4.03 X 10^^ photons ) 


At = 800 nm 


\p = 800 nm 


Ft = 1.2 


Fp = 4.5 


Wot ~ FtXt — 0.96/im 


"wop = FpXp = 3.6/im 


ZRt = TvwoJXt = 3.6/im 


ZRp = TTWop/Ap = 50.9/im 


Embedded physical phase based on QED 


S = 3.17 X 10"^ from ^ 




(Snq^d = 1.34 X 10"** from ^ with C = 7 and e = 7r/2) 




Focusing parameter in Signal Path 


expansion factor i?e = 5.0 x 10* 


Wp^ = ReWop = 18cm 


fp^ = 5m 


Peak intensity A^p^ = Ep/{2nw'^pJ ( 1.96 x lO" photons ) 



TABLE I: Parameters to produce Fig[5] and [6] based on the conceptual experimental setup in Fig|4l The subscripts t, p and 
Pe denotes parameters associated with the target laser, the probe laser and the probe laser after the expansion before the 
final focusing in the Signal Path in Fig|4j respectively. The definitions of parameters are explained in the text in section [Till 
and subsection IIIIBI As a basic constraint on the common beam diameter d of target and probe lasers in FigU] we assumed 
d ~ 1 m. 



several-photon sensitive photo device with an electrical 
amplification process such as the two-dimensional image 
intensifier camera. This device can increase the sensi- 
tivity to the number of countable photons and results 
in a higher resolution on the phase measurement. The 
peripheral sampling within 1cm is possible by locating 
flexible fiber bundles much like endoscopes with the ^m 
precision, as alr eady performed in a test experiment for 
another purpose [41 1. 



Integrated number of photons per pixel along y-axis per single shot 



C. Background in the phase contrast Fourier 
imaging 

The dominant background source of the current mea- 
surement may be caused by the refractive index shift due 
to the plasma creation from the residual gas along the 
path of the focused target laser pulse. The refractive in- 
dex of the static plasma in the limit of negligible collisions 
between charged particles is expressed as 




N 



1- 



(30) 



where ujq is the angular frequency of the target laser, ujp 
is the plasma angular frequency defined as y^ATre^ne/rrie 



and 7 is the relativistic Lorentz factor given as \/l + 
with ao = 0.85 x 10-^X[fim] ^h\Wlcm\ In the low 
pressure limit of the residual gas the amount of refrac- 
tive index shift A7V is expressed as uj^ j^'^uj^' . Although 
the refractive index in the plasma becomes smaller than 
that of the peripheral area with neutral atoms, the in- 
verted phase contrast of the phase shift inside the probe 



FIG. 6: The integrated number of photons over each 50 x 
50/xm^ pixel area along y (top) and i-axis (bottom) at the 
focal plane per single shot focusing. These correspond to the 
expected photon yields in a single pixel column along the y 
and X-axis in Fig[S] The plotted region is further extended 
to ±1 cm in both axises. All relevant laser parameters are 
summarized in TabU The ideal Gaussian pedestal distribu- 
tions without physical phase shift are also superimposed on 
the same figures. The number of photons at the focal spot 
are truncated to zoom the photon yield due to phase shift. 
By sampling photons in regions sufficiently distant from the 
focal spot, one can increase the signal-to-pedestal ratio. 
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pulse still maintains the rectangular shape along the tra- 
jectory of the target laser. Therefore, it should pro- 
duce the same characteristic diffraction pattern on the 
focal plane to that of the nonlinear QED case as ex- 
pected from the Babinet's principle. In order to re- 
duce this effect, we need to reduce the electron den- 
sity Ue in the residual gas. If we take 7 1 as the 
upper limit of AA^ estimate, the refractive index shift 
~ 10~^^ due to the nonlinear QED effect for a refer- 
ence energy density ^ IJ/zim"^, corresponding to the air 
pressure ~ 10~^ Pa. This pressure is easily attainable 
with conventional vacuum pumps. The collisional fre- 
quency due to interactions between electrons and ions is 
expected to be 10® — 10^s~^ at the critical electron density 
ncr[c-'rn~^] = 1.12 x 10^-'^/A^[^m] where cUp equals ujq. For 
duration time of ~fs of the target laser pulse, the inverse 
bremsstrahlung radiation due to the collisional process in 
the residual gas is negligible at ^ 1G~^ Pa. Therefore, the 
dominant background contribution from the residual gas 
plasma can be suppressed with the pressure well below 
10-6 Pa. 

We note that what actually happens is rather more 
dynamical due to the pondermotive force by the high 
intensity laser field. In such a case the refractive index 
shift based on the static plasma state gives only the upper 
bound on the amount of the local refractive shift. As 
discussed in IIII Bl we aim at the measurement of the 
refractive index shift on the order of ~ 10^® due to the 
nonlinear QED effect in the intensity assumed in TabU] 
This number is much greater than ~ 10~^^, as discussed 
above. Therefore, even if the dynamical effect is taken 
into account, it should be negligible to the measurement 
of the local refractive index shift, because the absolute 
value of the refractive index shift is small enough. 



IV. PROBING VACUUM FIELDS IN VERY 
LOW ENERGIES 

Extending on the QED investigation in the previous 
section, we now wish to explore the nonlinear vacuum 
response to lasers mediated by 'vacuum fields' whose en- 
ergies are much lower than the energy of the pump (and 
probe) lasers. To the current knowledge no such fields 
have ever been observed. This is in spite of the fact that 
there are numerous theories advanced. These theories in- 
clude: axions[23|, minicharged particles [i^ - lisf . and dark 
energy [2^ [231 . Why is this? This may be thought of as 
follows. After Rutherford's discovery of the inner core 
(i.e. nucleus) of an atom being very tiny compared with 
the size of the already tiny size of the atom, the exper- 
imental search went to explore ever smaller constituents 
of matter and thus the thrust for higher energy or mo- 
mentum experiments. Theories have gone hand- in-hand 
with this exploration, succeeding in ever shorter ranged 
interaction theories and unification of forces, typified by 
the electroweak theory [i^. We have referred to this 
standard and extremely successful method as the high 



momentum approach. Almost all laboratory research ef- 
forts have been on this approach to date. This approach, 
though successful in exploring high energy physics, is not 
suitable to explore energies much lower than eV. As men- 
tioned in SeclH these fields that might exist in much lower 
than eV cannot strongly couple to matter, because if it 
did, it would have been long discovered in lower ener- 
gies. Thus these fields, if they ever exist, must couple 
weakly. This means that we need to have an extremely 
strong driver to manifest a sufficient signal overcoming 
this weakly coupling interaction to show up above noise. 
We have not possessed sufficiently powerful such photon 
sources to date. Perhaps this may be changing now, how- 
ever, with more intense laser tools are to become avail- 
able fd\ . What we have called the high amplitude or high 
field method [l^, |Tl[ may provide an alternative path to 
detect such low energy weakly interacting fields that are 
spread over semi-macroscopic scales. 

When the constituent's energy is much lower than the 
probing photon's, it is our suggestion that we employ 
two laser beams in the co-propagating geometry. As we 
have motivated our method in the introduction, two co- 
parallel beams produce a very low center of mass energy 
interaction at the beat frequency (equal to the subtrac- 
tion of the two laser frequencies) [42]. The approach is 
a standard Brillouin scattering of laser in matter [1- 
m, [5^. In media such as condensed matter and 
plasma the utility of these two points have been recog- 
nized well in the study of nonlinear optics. For example, 
the Brillouin forward scattering (BFS) method [sj relies 
on two co-propagating laser beams with frequencies very 
close to each other, where the difference between these 
two is matched with the eigenfrcqucncy of the acoustic 
modes of the medium, typically one of the lowest eigen- 
modes. (We note that BFS is out of a more broad class 
of physical processes of the parametric instabilities [Slj . 
In addition a class of interaction under the modulational 
instabilities [57| shares a similar feature we may need.) 
It was known that this process allows strong coupling 
of photons to low frequency eigenmodes of the medium. 
This interaction could resonate with the constituent's 
very low eigenfrequency, should there be an eigenmode in 
its vicinity. Second of all, the co-propagating setup allows 
us to make the two beams interact over much prolonged 
interaction time, thus much amplifying the nonlinearities 
and signal arising from these. 

In order to pick up the experimental signal of strong 
coupling to the long-range mode, however, we suggest us- 
ing the second harmonic generation unlike Brillouin scat- 
tering where the detection is normally performed via the 
first order phonon-photon coupling. The pioneering re- 
search by Franken et al.[l| detected the nonlinearity in a 
quartz crystal. Two photons in co-propagation strongly 
interacted through the quartz fields over the coherence 
volume to produce the second harmonic generation. This 
process may be schematically looked upon as the case 
in FigjT^a). There the quartz nonlinearities have mixed 
two forward propagating photons (w) to produce a pho- 
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FIG. 7: Schematic diagrams of photon-photon interactions 
in matter and in vacuum, (a) Second harmonic generation in 
the experiment by Franken et al. l], arising from the non- 
linearity of crystal fields irradiated by (intense enough) laser 
fields; (b) Probing the QED vacuum nonlinearities, as sug- 
gested by Heisenberg and Euler [3ll . |3^ . where the vertices 
of the coupling in the Feynman diagram are characterized 
by the fine structure constant of vacuum a and thus a weak 
nonlinearity requiring us to employ much more intense fields 
than the case (a). The leading order interaction is the elas- 
tic photon-photon scattering, though as a higher order there 
exists a second harmonic generation as well; (c) Probing po- 
tential light-mass m fields in vacuum with intense laser fields. 
The vertices are characterized by very feeble couplings of M"^ 
and g [23. The expected second harmonic generation may be 
said to be not dissimilar to case (a). In order to increase the 
observable signal, suggestions have been made. 

ton with 2uj (and possibly another photon with frequency 
'--^ 0) . In the previous sections we have considered nonlin- 
ear QED process (FiglTl^b)). In this two incoming pho- 
tons are mediated by virtual electron-positron fields and 
outgoing are two photons. The extreme forward scat- 
tering amplitude with quasi-parallel incident photons is 
known to be largely suppressed in the QED process as 
we discuss later. This is because the center-of-niass en- 
ergy of the colliding two photons is too low to satisfy the 
relevant mass scale of electron-positron pair. 

In this section we consider a kinematically similar, but 
much more daring measurement of the photon-photon 
interaction mediated by low frequency "mode" of vac- 
uum by setting up intense co-propagating lasers. As we 
discuss later in detail, the process we focus is based on 
FigEl^c) which shows the second order in the photon- 
"mode" coupling. It is as if we are extending forward 
Brillouin scattering with the strong coupling to low fre- 
quency "mode" with the help of the resonance feature 
and Franken's method [l[ to obtain the clearer signature 
to probe the second order photon-" mode" coupling. 

In the previous sections, we have discussed the photon- 
photon interaction based on the QED process projecting 
for a first experimental verification of the real photon- 



photon interaction in the optical wavelength range. In 
addition to the verification, the ratio between the first 
and the second term in the one-loop effective Lagrangian 
in ^ can yield a general test to see whether vacuum 
contains other effects beyond QED examining its value 
at 4:7. Scalar field (p and pseudoscalar field a may con- 
tribute to the first and second term, respectively. Light 
scalar fields as a candidate of dark energy have been re- 
cently intensively discussed [26[, while the pseudoscalar 
fields (axion-like-particles) may be a source of dark mat- 
ter [i^l and also possibly dark energy [25|. However, 
we recognize that the test in the phase contrast Fourier 
imaging is very limited in the mass and coupling of those 
fields. In this section, therefore, we extend our method 
to search for those new types of fields by instituting co- 
propagating laser beams. These endeavors may be looked 
upon as in FiglT^c). Again two parallel photons come in, 
while two parallel photons come out. We note that our 
approach is similar to but distinct from many laboratory 
experiments with lasers (sj, Issl - lssj already performed 
and proposed to search for those light fields. 

Consider more details of the effective interaction La- 
grangian as illustrated in the triangle part in FiglTJ^c), 
where and cr couple to the electromagnetic field via 
quantum anomaly-type couplings [l^l 

-L„=g„M-^^F^,F^''a. (31) 

Here gM~-^ denoted by subscript (p ^i^d a for scalar 
and pseudoscalar fields, respectively provides the cou- 
pling strength. The dimensionless coupling g is typically 
proportional to the dimensionless fine structure constant 
a for the two photons to couple to the virtual charged 
particle pair in the triangle part. The effective coupling 
includes the large mass scale M to couple to the light field 
with the mass m. The large M induces the weakness of 
the coupling via M~^. For example, the Newtonian con- 
stant G is expressed as SttG — hcA-Ip'^, where Mp is the 
Planckian mass of 10^^ eV. The weakness of G is the 
manifestation of the large mass scale at the vertex in the 
triangle coupling. In what follows we abbreviate the sub- 
scripts (j) and a on the coupling gM~^ and the mass of 
light field m unless we need to explicitly distinguish the 
type of the fields. We use natural unit h = c — I through- 
out subsequent sections, unless we explicitly note. 

As a quite challenging case we have attempted a the- 
oretical approach to search for an extremely light scalar 
field via the resonance in [64?] , which may be sensitive to 
the mass scale of ~ 10~^ eV with M^^ = Mp^ cor- 
responding to the gravitational coupling. The method 
provides a new window into scoping the physics in the 
Planckian mass scale by photon interactions in the quasi- 
parallel incident laser beams in laboratories. In this sec- 
tion we review the essence of the approach and further 
develop basic formulae to apply the method to a concrete 
experimental setup in order to discuss reachable limits on 
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the coupling strength gM~^ and the mass m for a given 
laser intensity in the next section. 

As illustrated in FiglHl we introduce an unconventional 
coordinate frame in which two photons labeled by 1 and 
2 sharing the same frequency are incident nearly paral- 
lel to each other, making an angle § with the common 
central line along the z axis. We define the z — x plane 
formed by pi and p2- The components of the 4- momenta 
of the photons are given by pi = {ujsm'd,0,ujcos'd;uj) 
and the same for p2 but with the sign of -i^ reversed, 
and P3 — (cj3 sin 63,0,003 cos ^3; W3) and P4 with UI3, 03 re- 
placed by CJ4, —64, respectively. The angles 63 and 64 are 
defined also as shown in Fig|8l This coordinate system 
can be transformed to the center-of-mass (CM) system 
for the head-on collision {•& = tt/2) by the Lorentz trans- 
formation with u/c — > 1 for i? — ;> 0. Conversely, this 
implies that the realization of the quasi-parallel collision 
in the laboratory frame corresponds to the realization of 
the extremely low CM energy, as we see below. 

In this frame one of the final photons in the forward di- 
rection along the z axis must have an upshifted frequency 
due to the energy-momentum conservation independent 
of the physical origin of the dynamics. In the limit of 
1? — 0, a process of 103 — 2a; is realized. This frequency 
doubling nature is an extremely valuable characteristics 
from the experimental point of view, as compared to the 
case with no frequency shift in the centcr-of-mass system. 
In addition, more importantly, it is essential to maintain 
a quasi-parallel nature of the incident beams in order to 
access the resonance point, as we shall stress later. 

The energy-momentum conservation laws are 

0-axis : 103+0:4 = 2oj, (32) 
z-axis : a;3 cos 63 + 004 cos 64 = 2oj cos (33) 
a;-axis : W3sin^3 = aj4sin(?4. (34) 

From the conditions < a;3_4 < 2oj, we may choose < 
63 < d < 64 < TT, without loss of generality. From ((32)) - 
we derive the relation 



sm 6^3 = sm 04 



sin^z? 



1 — 2 cos 1} cos 04 + cos^ 1) 



(35) 



The differential elastic scattering cross section favoring 
the higher photon energy a;3 is given by 



da 



-^sm-^d{oj3/2ojf\M\'', 



where Ai is the invariant amplitude and 

oj sin^ i9 



0J3 



1 — cos z9 cos 63 



(36) 

(37) 
2oj, as 



Here we expect the upshifted frequency 0J3 
(?3 ^ for 1? — ;> 0, as mentioned before. 

The resonance decay rate of light fields with the mass 
m into two photons is expressed as 



P2 e(l 




FIG. 8: Definitions of kinematical variables for the suggested 
co-propagating photons (this figure is quoted from fi4^). 



The light field is exchanged between the pairs {pi,p2) 
and {p3,P4), thus giving the squared four- momentum of 
the scalar field 

g2 ^ (p^ + p^f = 2cj2 (cos 2d - 1) (39) 

with the metric convention (+ + H — ) for the definition 
of four momenta. 

With the polarization vectors given by e^*-'^'' where i = 
1, • • • , 4 are the photon labels, whereas /3 = 1, 2 are for 
the kind of linear polarization as depicted in Fig|Sl we 
summarize the non-zero invariant amplitudes for scalar 
exchanges 

Mini = A^2222 = -X1122 = -A^22ii, (40) 
and for pseudoscalar exchanges 

A^i2i2 = A^i22i = -X2112 = -A^2i2i, (41) 

where the first two digits in the subscripts correspond 
to the states of the linear polarization of incoming two 
photons 1 and 2, respectively and the last two correspond 
to those of outgoing two photons 3 and 4, respectively as 
iUustrated in Fig El 

We focus on one of these non-zero amplitudes by de- 
noting it as A4; 



M 



-{gM- 



-1\2 



(cos2z9- 1)^ 
2cj2 (cos2i?- 1) -|-m2' 



(42) 



where the denominator, denoted by D, is the light field 
propagator. We note qs is timelike. We then make the 
replacement 



(to — iF) K, — 2imV . 



(43) 



Substituting this into the denominator in (1421) . and ex- 
panding around to, we obtain 



(38) 



V 



-2 (1 - cos 2-3) (x -I- m) , with x = - w^, (44) 
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where 



toV2 
1 - cos 2-3 ' 



1 - cos 2-d 



From (pS)) and pS)) . a is also expressed as 



8"7r 



(45) 



(46) 



which exphcitly shows the proportionaUty to M^"^ . We 
then finally obtain the expression for the squared ampli- 
tude as 



\M\ 



(47) 



As for the off-resonance case % 3> a, |A^p is largely sup- 
pressed due to the factor a? oc for the case of small 
coupling M"-"^. On the other hand, if experiments take 
ideally the limit of w — s> w^, — > (27r)^ is realized 
from (|T7)) . This is independent of the smallness of the 
factor M^^ as expected from the off-resonance case or 
equivalently from the square of (|42|. This is the most 
important feature arising from the resonance that over- 
comes the weak coupling stemming from the large rele- 
vant mass scale such as M — Mp. However, we then con- 
front the extremely narrow width a for e.g. gm <^ 1 eV 
, M Mp = 10^'^ eV and Wr ~ 1 eV. We now discuss 
how to overcome this difficulty. 

In conventional high energy collisions the beam mo- 
menta are implicitly supposed to be their mean values. 
This is because the momentum spread, or the uncer- 
tainty expected from the de Broglie wavelength of the 
relativistic particle is negligibly small compared to the 
relevant momentum exchanges in the interaction that ex- 
periments are interested in. Resonance searches in such 
experiments must adjust x in (011) and (|47p so that the 
mean x is close enough to the peak location within ±a. 
On the other hand, in the case of the co-propagating laser 
beam aiming at detection of extremely small momentum 
exchanges via the resonance, the situation is quite differ- 
ent due to the nature of incident waves. This is because 
the uncertainty included in the initial photon momenta 
is much larger than the relevant energy scale of the reso- 
nance, leading to the condition |x| a. In this case the 
squared scattering amplitude must be integrated over the 
possible uncertainties on the incident wave function on 
the event-by-event basis. As we will explain below, the 
uncertainty in the incident photon momenta is related 
to the uncertainty in x via the uncertainty in the inci- 
dent angle d. Therefore, it is instructive to consider the 
feature of the integral of the resonance function of the 
Breit-Wigner(BW) formula j66|] from x- to x+ as fol- 
lows; 



Tdx 



a tan ^— ^ 



X- 



(48) 



where / = a7r/2 and I = air for x+ — ~X- — and 
X+ = —X- — oo, respectively. This indicates that the 



value of integral is proportional to a, i.e. from (|46p . 

The value ranges for the finite and infinite integrals over 
only factor of two. From this fact, we expect that the 
integral enhances the squared scattering amplitude by a 
factor of compared to the non-resonance interaction 
proportional to A/^* from (I42|) . as long as the peak is 
contained within the experimental resolution on x, i.e. 
the condition x+ > a and X- < ~o is satisfied. This im- 
plies that experiments in the case of the co-propagating 
laser beam configuration do not have to make efforts to 
adjust X close to the extremely narrow resonance region 
thanks to the consequent huge enhancement by the inte- 
gral over the wide range on x- Meanwhile, it is difficult to 
identify the exact location of the resonance mass within 
the wide gate on x- Since we are interested in having the 
sensitivity to the extremely weak coupling such as grav- 
ity M = il/p, the enhancement of the squared amplitude 
is more crucial than finding the exact location of reso- 
nance masses. As we discuss in the following sections, 
however, we may be able to provide a crude estimate on 
the order of the mass scale of the resonance even in such 
a situation. 

The consideration above leads us to parametrize the 
squared scattering amplitude as follows: 



4-00 



p{pi,P2)\M\^dx 



(27r) 



+ 00 



PiPl,P2) 



X^{pi,P2) +a' 



:dx, (49) 



where p{pi,p2) is the normalized probability distribution 
to supply nominal combinations {pi,P2) from incident 
laser fields and the real part x(pijP2) is indirectly speci- 
fied by {pi,P2) via the incident angle between incident 
two photons. This parametrization expresses averaging 
over possible combinations of pi and p2. We note that 
(pi J P2) is not a priori observed momenta, but just a nom- 
inal specification among the possible initial momenta. In 
other words, p{pi,p2) is not a statistical weight on the 
discrete momenta after the contraction of the wavepacket 
of each photon state. This implies that an infinite statis- 
tics is not necessary to obtain the continuous nature. 
Rather, we need this treatment even for a two photon 
state as long as the source of photons is not the perfect 
plane wave. This allows the continuous integral on x via 
the continuous combination of p{pi,P2) in (|49p . As long 
as the probability weight p(pi , P2) is close enough to unity 
around the resonance peak, the enhancement discussed 
with is guaranteed. This is the essence of our main 
strategy of the co-propagating configuration to overcome 
the difficulty due to the narrow resonance width a. 

In order to design experiments, we start from the res- 
onance condition, the first of (j45p . by assuming z9 ^ 1, 



2i!)u 



(50) 



We note that the product 2'duj corresponds to the CM 
energy of incident two photons. This indicates that ex- 
periments have the two adjustable handles for a given 
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mass scale or the CM energy. We emphasize that we can 
lower the CM energy by several orders of magnitude by 
only introducing smaller with fixed lo. This should be 
contrasted to high energy colliders where we need a large 
amount of efforts to increase the CM energy by an order 
of magnitude. This advantage is also supported from a 
technical point of view, since scanning the incident angle 
d should be much easier than scanning the energy uj of 
the resonance. We point out that the resonance condition 
([50]) is not just at one point, but rather in a hyperbolic 
band in the i) — ui plane given with a finite resolution 
with 5d of the incident angle -i^. This implies that the 
deviation 5uj from the resonance energy can satisfy 
the same resonance condition with a different § within 
■iiSd. As far as Soj/u) <^ S-d/d is satisfied in the setup, we 
can ignore the effect of Slo. This is in fact the case, as 
seen in the following discussions. Therefore, we can take 
the attitude that we fix the incident energy at the optical 



frequency uj 



■'opt 



and scan m by changing d around dr , 



where ujopt and i?r satisfy the resonance condition based 
on (EOl 



= mV(4i?2). 



(51) 



In this case the variation on § leads to the variation on 
X- This is expressed by the following relation based on 
the second of 



'■^opt 



4^2 ^opt 



(52) 



where £ = d /"dr in the unit of dr is introduced. 

We now discuss the average of the squared amplitude 
jyVlp over the possible uncertainty on the incident angle 



-it/2 







7r/(2iJr) 



p{e)\M\^drde, 



(53) 



where p{'d) is a probability distribution function normal- 
ized between and 7r/2 as a function of continuous un- 
certainty on 'd from arbitrarily chosen two photon com- 
binations within a laser pulse. The incident angle is 
re-expressed with e with dd — "drde for the second of 
([55)1 . From the relation ([SI]), we find e = (1 -x)"^/^ ^-^^^ 
de = 1/ [2u)1pf)e^dx. Equation ([55)) is then re-expressed 
with X 



|7U|2 = (27r) 



dx. (54) 



where ([17|) is substituted. This equation is the exact 
representation of (I49|) starting from the uncertainty on 
the incident angle t9, if the upper limit of the integral 



range is regarded as large enough compared to a. Let us 
define a: = to explicitly discuss the structure of the 
integral kernel in the unit of the width a of BW. With ^, 
([M]) is further re-expressed as 



|7W|2 = {2Tiy 



(1 - 0^3/2 ^2 ^ I 



dt (55) 



where the first factor of the integral kernel corresponds to 
a normalized weight function and the second is BW with 
the width of unity. This expression explicitly shows the 
enhancement by the factor of a implying the proportion- 
ality to M~2 based on ([46|l . As long as p is a monotonic 
function, the weight function in front of BW can be close 
to unity for small ^ because of ^ 1. With such a 
weight we expect that the value of the integral may be 
close to that of the BW as we discuss with a concrete 
weight function in the following section. 

Remaining issue is how to further cope with the prob- 
lem of still very small M~^, though much larger than 
A/^"*, in experiments. First, this can be solved by sin~^ 
behavior of the cross section in ([55]) that arises from the 
phase volume factor and the flux factor in the quasi- 
parallel two photon interaction. For extremely light 
mass, this factor gains a large number due to small 7?^- 
Second, the intense laser fields can provide a large lumi- 
nosity and the intensity of the signal is proportional to 
square of the intensity of laser in the limited case of the 
incoherent two photon interaction. We have three ingre- 
dients or knobs; the M'^ enhancement by the weighted 
BW integral, the factor of and the growth of laser 
intensity. By marshalling these knobs, we expect to in- 
crease the detectability for the undiscovered light fields 
in vacuum, which have evaded from our grasp to date. 

In the following sections we consider the experimen- 
tal realizations with ujr ^ 1 eV (optical laser) and 
0(1) < i^r < 0(10^^*^) by scanning of the mass range 
0(1) < m < 0(10^1°) eV. We then plug explicit weight 
functions into (|55p based on the suggested experimental 
setup. By combining and ([M]) . we discuss reachable 
mass-coupling limits for a given laser intensity attainable 
in future experiments. 



V. CONCEPTUAL LASER MEASUREMENT 

We emphasized the importance of photon-photon in- 
teraction in a co-parallel or at small angle setup in or- 
der to increase the signal of detecting low energy con- 
stituents. Thus a simple way is to explore the mass range 
m < TTUJ by using two independent laser beams with the 
small incident angle. We then directly measure the res- 
onance curve in (j50[) by scanning both and uj to quan- 
titatively observe the nature of the resonance curve. For 
the very smaller mass scale, or equivalently smaller inci- 
dent angle, however we must take into account the beam 
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spread in the diffraction limit. This determines the con- 
trohable smallest incident angle, or the mass scales of the 
light fields we look for. We consider here one-beam fo- 
cusing in order to provide the simplest basis to quantify 
reachable mass-coupling limits for given a set of experi- 
mental parameters. 

The conceptual experimental setup with one-beam fo- 
cusing geometry is illustrated in FiglQl Incident photons 
from a Gaussian laser pulse with linearly polarization 
are focused by the conceptual lens component into the 
diffraction limit. Quasi-parallel incident photons inter- 
act with each other between the lens and the focal point, 
from which photons 3 and 4 are emitted nearly in the 
opposite direction along the z axis with W3 ^ luj and 
a;4 ^ 0. The mirror with dichroic nature is transpar- 
ent to the non-interacting photons with the beam energy 
of while ^3 is reflected to the prism (equivalent to a 
group of dichroic mirrors) which selects ^3 among resid- 
ual io and sends it to the photon detector placed off the 
z-axis. This process is assisted by the polarization fil- 
ter. From the polarization dependence of the invariant 
amplitude in (^(1)) and (HI]), the combinations of polariza- 
tions of two photons between the initial and final states 
must satisfy 11 — >■ 11(22) for the scalar field exchange 
and 12 12(21) for the pseudoscalar field exchange. We 
note that we can choose the type of fields we search for 
by setting the initial polarization state. In the case of 
one-beam focusing the search for the scalar field is eas- 
ier, because we do not have to mix the two polarization 
states as in the case of the pseudoscalar field. Further- 
more, the selection of the rotated final state 22 can en- 
hance the signal-to-background ratio for the scalar field 
case, because a huge number of non-interacting photons 
has the final polarization state of 11. In what follows we 
provide formulae to evaluate the accessible limit on the 
mass-coupling defined in pip for a given laser intensity, 
in the case that we detect a double frequency photon per 
laser focusing. 

The Gaussian laser parameters are summarized in Q 
through (|13l) . We now estimate the effective luminosity 
L over the propagation volume of the laser pulse. We 
now restore physical dimensions of 1% and c in this section, 
unless we explicitly note. Consider a Gaussian laser pulse 
with duration time r with the speed of light c and the 
average number of photons N per pulse. The exchange 
of light field may take place anywhere within the volume 
defined by the transverse area of the Gaussian laser times 
the focal length / before reaching the focal point. We 
first consider the effective number of photons iVi„t during 
an interaction with the time scale of At. As a result of 
the interaction we observe a frequency doubled photon in 
the laboratory frame. The momentum transfer of ^ hio j c 
between photons defines the minimum interaction time 
scale from the uncertainty principle as follows 



as 



(56) 



The effective number of photons during At is expressed 



At 



-N. 



(57) 



Making the pulse duration r ^ At maximizes the instan- 
taneous luminosity. Suppose a point z along the laser 
propagation axis. The instantaneous luminosity at the 
point z is defined as 



L{z) 



Hz) 



27rWn ^2 + Zr 



(58) 



where C{Nint,'2) denotes a combinatorics to choose two 
photons among a large number of photons available 
within time scale of At and the expression of w'^{z) in 
is substituted to obtain the second with the approx- 
imation for the combinatorics. We then consider the av- 
eraged instantaneous luminosity L over the focal length 
/ as follows; 

L = [\{z)dz - 

-Jfl^ZRtan-\f/zn) = -f^tan-i(//zfl,). (59) 
2Trfw^ 2/A 

The number of effective bunches b is related with / as 
. / 



cAt' 



(60) 



The effective luminosity C over the propagation volume 
of the laser pulse is finally expressed as 



C^bL 



At 

~2ctX 



tan-i(//zfl),(61) 



where (1571) is substituted in the last step. 

The uncertainty on the incident angle between two 
light waves is expected to be 



Ad 



zr 



(62) 



from the definition of the Rayleigh length zr = ttwq/X. 
The minimum beam waist wq at z = varies with the ex- 
perimental conditions, the focal length /, Rayleigh length 
Zr and diameter of conceptual lens d, 



Wo = id/2) 



if/zR) 



V^ + if/zny 



(63) 



We see that Ad may be controlled via wq by choosing 
suitable / and d in experiments. 

The possible uncertainty on the inci dent an gle i9 affects 
the average of the squared amplitude as shown in 

the first of (|53p . In order to obtain an approximation 
close enough to reality, we plug the following step func- 
tion into (l53l): 




for < i9 < At? 
for Ad <i9 < n/2 



(64) 
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yield per laser pulse focusing as follows: 



P2 




polarization 
7 ' filter 

photon detector\3'* 



FIG. 9: Suggested experimental setup of the co-propagating 
photon interaction and detection. The linear polarizations of 
incident and outgoing photons are drawn only for the scalar 
exchange with the scattering amplitude |A^ii22| as an exam- 
ple. 



which is normalized to the physically possible range < 
I? < tt/2. By substituting (|M)) into we obtain 



_ (27r)^ l^r 
a-i{l-(i?,,/Ai5)2} 



(1 - aO^/^ (e + 1) 



de, (65) 



where the first factor of the integral kernel corresponds to 
the weight function and the second is Brcit-Wigner(BW) 
with the width of unity. The weight function of the ker- 
nel is close to unity for small ^ due the smallness of a 
in (|46l). Therefore, the value of integral in ([65)) is al- 
most equivalent to that of B W [g^ . This is because the 
monotonic positive weight function approaches zero as 
^ — >■ — cxD more rapidly than the pure BW, whereas the 
pure BW suppresses the increase of the weight function 
close to zero at ,^ ^ a-^{l - {dr/A'df} for Ai? < 1. We 
then approximate (p5)) as the integrated BW over ± ~ 00 
as follows: 



(27r)^ 



(66) 



With a in (|46|) and in ([66]) . the differential cross 

section in p6p is finally expressed as 



/ da 



TT /27r 

64 It 



Is) (f 



where the approximations UJ3 2Tiijj and ??r ^ 1 are 
taken into account. 

Multiplying (|5T|) by (|57)) . we obtain the differential 



dy _ f da 

dns ydfia 

1 At X 
5127r r ct 



tan ^if/zn) 



Ai9 \M J 



There are several experimental knobs to affect the ob- 
servable events in ([68]) . If we choose t ^ At ^ X/c 
resulting in cr '-^ A, we can maximize the effective lumi- 
nosity. From (|62l) , the reduction of A-d or increasing wq 
enhance the yield. From ([55]). we express f/zn as 



{d/2y 



(69) 



From this relation, increasing wq as large as (c?/2) also 
enlarges f / zr. This introduces a slight increase for the 
factor tan~^(//zH), though its effect is tempered by the 
nature of tan^^. 

As a short summary, we make the most important note 
from the experimental point of view based on this con- 
ceptual design. The condition dr/ Ad — 1 maximizes the 
chance to search for the resonance, while dr/ Ad > 1 re- 
sults in a huge suppression of the cross section by M""* 
(h. — c = 1) as we discussed. This is because the res- 
onance peak is out of the region covered by Ad. This 
parameter corresponds to the sharp cut-off of the cross 
section. Therefore, controlling Ad via relations (|62|) . ([63)) 
and (|69p can provide an experimental way to define the 
mass range that we eliminate if no signal is found. From 
the resonance condition in (|50p . we evaluate the lower 
bound on the mass range wc can exclude by this concep- 
tual design as follows (?i = c = 1): 

rrimin ^ 2drLUopt ~ 2AduJopt 

^ 2Tr-\X/wofuJopt > 2Tr-\X/dfiJopt, (70) 



where (|63[) is substituted and the obvious experimental 
condition wq < dis required for the last inequality. Below 
this lower mass bound we cannot discuss what the mass 
scales of the light fields are, even though the measure- 
ment is still sensitive to the presence of lighter resonances 
than the lower mass bound due to the large enhancement 
by d-^. 



VI. DISCUSSION BASED ON ONE-BEAM 
FOCUSING 

As a demonstration we now discuss the necessary laser 
intensity for the following reference case based on the 
one-beam focusing setup. First, we note that At in ([55]) 
is the resolvable minimum time scale. As long as we dis- 
cuss extremely low mass field, the interaction time scale 
may be over h/mc^. In the case of low mass the inter- 
action time scale may be much longer than 27raj^(. On 
the other hand, photon-photon interactions must occur 
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within r. This fact imphes we should assume At/r = 1 in 
the case At > r. Since the effective luminosity tends to 
be enhanced for smaller beam waists, we assume a short 
pulse laser with duration time close to 2Truj'^p^ in the fol- 
lowing discussion. Let us assume that the diameter of the 
conceptual lens is ~ 2 m by taking the damage threshold 
of the lens into account for the use of high intense laser 
pulse. With the assumption / <C zr, the minimum beam 
waist in ([63]) is approximated as wq ^ ((i/A/(27r))^/'^. 
This gives the cut-off of the sensitive mass region mcut 
defined by the condition < A-d via the relation (p^ 
as follows; 

^ 2u;„,,A^ = ^"°^ ^(d/A/(27r))^/3 ' ^^^^ 

If we are allowed to change the focal length between 
~ 1 < / <~ 1000 m, the cut-off mass range varies from 
~ 10^^ to ^ 10^^^ eV with ujopt ~ 1 eV. By choosing 
the focal length realizable in laboratories as / 3 m, we 
expect rricut ~ lO"^'^ eV. For uJopt ~ 1 eV, the resonant 
incident angle is -Or ^ 10^^° from ([50]) . As the most chal- 
lenging case, we assume the coupling as weak as gravity 
Mp^. We are now ready to estimate the aver- 
age number of photons iVi to expect dy/dU^ ~ 1 per 
pulse focusing based on (1551) for the physical parameters: 
TO - 10-1" eV, g ~ a = 1/137 and Mp - 10^^ eV, 
and for the experimental parameters discussed above: 
ijJopt ~ 1 eV, r ~ 10 fs, At/r = 1, wq ^ 0.01 m, d ~ 2 m, 
/ - 3 m, A?9 - 4 X 10"^ > i^r and f/zR - wq from §^ 
resulting in tan^^{f / zji) wq. By taking all the factors 
into account, we find Ni ~ 0(10^^), corresponding to 
^ 10 kJ per pulse focusing. 

If the one-beam focusing setup cannot detect any sig- 
nals, it excludes all of mass range below rricut and all of 
couplings stronger than gM~^ in that mass range. This 
is because, in principle, the sin""* "dr dependence in ([36]) 
enhances the cross section for the entire lower mass range 
below rricut- We note that there is no known physical 
cut-off on the minimum § in the photon-photon inter- 
action. This by itself can be a subject to be studied, 
since it is possibly related with the texture of vacuum 
which may prevent the smooth propagation of photons 
in vacuum. This may introduce a finite cut-off in the 
minimum incident angle between two photons. If the in- 
tensity discussed above is available, we may be able to 
search for extremely light fields in the totally unprobed 
domain: to < 10"!° eV and gM"! > lO'^^ GeV-^. 
This is a remarkable improvement compared to the ax- 
ion (pseudoscalar field) searches so far taken. They pro- 
vided the upper limit in the domain: m > 10~^ eV 
and gM-^ > IQ-^^ GeV'^ The way to estimate 

the mass-coupling limit in our suggestion should be con- 
trasted to the conventional ones 23] where experiments 
provided upper bounds on the mass-coupling limit. The 
gap between the two domains, however, can be filled by 
changing the experimental parameters in our concept. 
Even if there is no signals in the one-beam focusing, 
we only have to update the condition so that it satisfies 



Ai? > -dr for heavier masses by increasing Ad. In such 
heavier mass region, however, two-beam crossing geome- 
try relaxes the constraints on the optical design such as 
focal length. In either case the one-beam focusing setup 
considered in this paper provides a basis to define the 
mass-coupling limit as well as the necessary beam inten- 
sity as we have demonstrated here. 

Our arguments above have been based on the approach 
in which each interaction arise incoherently yielding the 
observed result by N^\M\'^. It is worth noting, how- 
ever, the interaction may occur coherently because the 
light field exchange gives a sufficiently long-range interac- 
tion, and hence individual scattering centers in the beam 
are indistinguishable. This coherence effect results in 
^|iVA^P which is enhanced by a factor of N over the 
yield expected from the incoherent summation, as was 
discussed in the long- wavelength neutrino interaction on 
a bulk target [68*1 . Upon more careful analysis of our sys- 
tem the question may be addressed if the required laser 
intensity can be relaxed due to a collective frequency shift 
rather than the shift on the single-photon basis. 

A major instrumental background for the frequency 
doubled photon is in principle expected to be the sec- 
ond harmonic generation (SHG) from the final focus- 
ing optical element. The dominant source of SHG may 
be the interface between the residual gas and the sur- 
face of the optical element where the centrosymmetry 
is maximally broken. Even from the maximal estimate 
~ 10^'^W/cm^ for a typical damage threshold, we ex- 
pect a negligible amount of 10~^° SHG photons from 
a Im^ aperture size with a 10 fs irradiation, if the op- 
tical components are housed in the vacuum containing 
10^° atoms/cm^ (~ 10"^ Pa) The confirmation of 

the negligible amount of the background SHG is a crucial 
subject for the present concept. 

As a dominant physical background we expect the 
lowest-order QED photon-photon scattering with the for- 
ward cross section ~ {a^ /m'^)'^uj^'d'^ [6^. This turns out 
to be much smaller than (p7|) due to the distinct behavior 
with respect to the incident angle d. This indicates the 
lowest-order QED contribution is negligible. 



VII. CONCLUSION 

We have suggested an approach to probe the nature 
of vacuum with intense lasers. Two main methods have 
been explored. The first is the phase contrast Fourier 
imaging at the focal plane to measure the phase shift of 
propagating light under the intense laser field. The sec- 
ond is the extremely light resonance search via higher 
harmonic generation by focusing a single intense laser. 
Both are based on similar ideas extended from those al- 
ready developed to probe matter. 

In these methods we take note of the nonlinearities 
of vacuum that are either considered to exist but at an 
extremely minute level or speculated to wait for our sen- 
sitive detection of an extraordinarily feeble signal. In or- 
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FIG. 10: Experimental ballparks of various approaches to 
probe matter and vacuum in the system size vs. the energy 
density. Selected systems are LC: electron-positron collision 
in the center of mass energy Ecms = ITeV at the future lin- 
ear collider [t^] assuming the electron size 10~^*cm, LHC: 
proton-proton collision in Ecms = 14TeV at Large Hadron 
Collider [zO], RHIC: gold-gold coUision in Ecms = 200GeV per 
nucleon pair at Relativistic Heavy Ion Collider [t^, the rest 
proton indicated by the asterisk as the origin of this plot, ELI: 
an optical laser pulse expected in ELI project [7l[, Schwinger: 
Schwinger limit [8|, Sun, the Milky Way Galaxy and event 
horizon with fltot ~ 1.0 and h = _H'o/100[km/s/Mpc]~ 
0.7 [T^I- The energy density axis is qualitatively interpreted as 
the inverse of the force range or the mass scale m of exchanged 
force, because the mean free path becomes shorter in higher 
density states as long as the coupling to matter is not weak. 
On the other hand the coupling strength to matter gM~^ de- 
fined in pi|l qualitatively refiects the necessary size of matter 
or vacuum in order to make the interaction manifest. The 
arrow to the higher energy density toward Schwinger limit is 
the direction to probe nonlinear QED interactions and also 
understanding of the non-perturbative nature of the intense 
field, while the arrows directing lower energy density region 
indicate the extensible domain by using co-propagating in- 
tense laser fields, since the mass range may be covered from 
~ 1 eV to ~ 0.1 neV and the coupling may be probed as 
weak as gravitational coupling for lighter mass scales. The 
energy density in this direction depends on the context. In 
the context of the scalar field as a candidate of dark energy 
in [2^ . Ig^ I , the energy density should be close to that of the 
event horizon. 



der to detect these weak nonlinearities, we need to spec- 
tacularly enhance the signal. The large leap in enhanc- 
ing these signals is garnered by the combination of (1) 
the rapid development of intense laser technology and its 
adoption here; (2) the employment of our suggested tech- 
niques that circumvent potentially huge noise and allow 
the enhanced interaction with the pursued fields. The 
former element (1) may be brought in, for example, by 
the intense optical laser beyond kJ. For the latter factor 
(2) we have suggested for the nonlinear QED problem a 
method of avoiding the pedestal noise in our phase con- 



trast Fourier imaging. For the exploration of possible 
new low mass fields we have suggested a method to hit a 
resonance with photons co-propagating over a long dis- 
tance. 

With the phase contrast Fourier imaging the nonlinear 
QED effects of the Euler-Heisenberg Lagrangian should 
be detected no more than the energy of laser around 10 J. 
The method provides a window for scoping vacuum via 
the dynamics of electron mass scale virtually. Such de- 
tection has never been made to date, and it heralds the 
research in the physics of vacuum with the high field ap- 
proach. With the detection of second harmonic genera- 
tion in the co-propagating setup we should be able to sur- 
vey a large sweep of the energy domain (eV to 10~^°eV) 
of the intermediating vacuum fields. If and when we pick 
up some signal in one particular energy range, perhaps 
we can zoom in to this specific energy (and thus wave- 
length) of photons by arranging the various knobs such as 
the crossing angle and the (long) beat wavelength of the 
electromagnetic waves by setting up a specific resonance 
cavity, we may be able to further increase the sensitivity 
and more deeply study their properties. 

Given the high intense optical laser beyond IkJ per fs- 
pulse duration in the near future, the realization of these 
suggestions may become an exciting challenge for future 
experiments of explorations of the physics of vacuum. 

Figure [TOl illustrates the experimental ballparks of var- 
ious approaches to probe matter and vacuum, in terms 
of the system size as a function of the energy density. 
The energy density axis is qualitatively interpreted as 
the inverse of the force range or the mass scale m of 
exchanged force in because the mean free path be- 
comes shorter in higher density state as long as the cou- 
pling to matter is relatively strong. On the other hand, 
the coupling to matter gM~^ in p2|) or ((3T|) qualitatively 
reflects the necessary size of matter or vacuum in order 
to make the interaction visible. The Galileo type tele- 
scope observes gravitational phenomena. These are on 
the extremely weak coupling of Mp^ with zero mass ex- 
change. High energy particle colliders, the Rutherford 
type microscopes focus on particle generation phenom- 
ena. These are due to strong couplings with heavy mass 
exchanges within fm scale. There is a huge gap between 
these two approaches. In other words the region of weak 
couplings with finite but light mass exchanges have been 
hardly probed so far. It is quite natural to start explor- 
ing if there exist important pieces of the puzzle of nature 
in these domains. These explorations might allows us 
deeper understanding of the nature of vacuum such as 
dark energy [73j . The progress of modern physics has 
been simply driven by those two types of the experimen- 
tal approaches. The proposed method with high intense 
laser probes semi-macroscopic vacuum compared to par- 
ticle physics and much smaller scale of vacuum compared 
to cosmology. Provided such semi-macroscopic vacuum 
scope, we increase our observational window into a new 
parameter regime of vacuum. 
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